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ABSTRACT 

We  apply  energy  methods  to  prove  an  existence  theorem  for  the  Dirichlet  initial 
boundary  value  problem  in  incompressible  nonlinear  elasticity. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  elastic  materials,  such  as  rubber,  can  be  subjected  to  large  volume-preserving 
deformations,  but  show  little  compressibility.  Such  materials  are  often  modelled  as  incom¬ 
pressible  elastic  solids.  In  this  paper,  the  authors  study  the  initial-boundary  value  problem 
for  the  differential  equations  describing  such  materials  and  prove  that  it  is  well-posed. 
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AN  EXISTENCE  THEOREM  FOR  THE  DIRICHLET 


INITIAL-BOUNDARY  VALUE  PROBLEM 
IN  INCOMPRESSIBLE  NONLINEAR  ELASTICITY 


William  J.  Hrusa1  and  Michael  Renardy2 


1.  Introduction  and  statement  of  results 

The  goal  of  this  paper  is  to  show  how  energy  methods  can  be  applied  to  establish  the 
local  (in  time)  existence  of  classical  solutions  to  the  Dirichlet  initial-boundary  value  prob¬ 
lem  in  incompressible  nonlinear  elasticity.  For  compressible  nonlinear  elasticity,  the  initial 
value  problem  posed  on  all  of  space  was  solved  by  Hughes,  Kato  and  Marsden  [7]  using 
methods  of  semigroup  theory,  and  the  Dirichlet  initial-boundary  value  problem  was  solved 
by  Kato  [8],  Chen  and  von  Wahl  [3]  and  Dafermos  and  Hrusa  [4],  Schochet  [14]  proves  an 
existence  theorem  for  the  initial  value  problem  in  the  incompressible  case  by  considering 
incompressible  materials  as  a  limit  of  compressible  materials  and  deriving  uniform  esti¬ 
mates  which  allow  passing  to  the  limit.  His  proofs  are  based  on  ideas  of  Klainerman  and 
Majda  [9,10],  who  considered  the  incompressible  limit  in  gas  dynamics.  A  direct  existence 
proof  for  incompressible  elasticity  (without  using  approximation  by  compressible  materi¬ 
als)  is  given  by  Ebin  and  Saxton  [5].  The  results  of  Renardy  [13]  can  also  be  applied  to 
incompressible  elasticity.  In  this  way,  an  existence  theorem  for  the  incompressible  initial 
value  problem  is  obtained,  but  it  is  not  proved  that  solutions  of  the  incompressible  problem 
are  limits  of  solutions  for  the  compressible  case.  To  our  knowledge,  there  are  no  previous 
results  on  initial-boundary  value  problems  for  incompressible  elasticity. 

In  all  the  papers  quoted  above,  energy  estimates  play  an  essential  role.  A  first  order 
energy  estimate  can  always  be  obtained  by  considering  the  physical  energy  of  the  elastic 
body,  which  is  the  integral  of  an  expression  involving  first  derivatives  of  the  unknown  de¬ 
formation  function.  Bounds  on  this  energy,  however,  are  not  strong  enough  to  be  useful  in 
an  existence  proof.  The  idea  is  to  differentiate  the  equations  of  motion  and  consider  the 
analogues  of  the  energy  for  higher  derivatives  of  the  solution.  The  a  priori  bounds  thus 
obtained  are  then  used  in  a  contraction  argument  establishing  the  existence  of  a  solution. 
There  is  an  essential  difference  here  between  the  initial  value  problem  on  all  of  space  and 

1  Department  of  Mathematics,  Carnegie-Mellon  University.  Pittsburgh,  PA  15213. 

2  Department  of  Mathematics  and  Mathematics  Research  Center.  University  of  Wis¬ 
consin.  Madison,  W1  53706.  After  August  1986:  Department  of  Mathematics,  Virginia 
Polytechnic  Institute  and  State  University,  Blacksburg,  VA  24061. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DA AG29-80-C-0041,  by 
the  National  Science  Foundation  under  Grant  No.  DMS-8451761  and  by  the  United  States 
Air  Force  under  Grant  No.  AFOSR-85-0307. 


the  initial-boundary  value  problem.  For  the  initial-value  problem  one  may  differentiate 
the  equations  with  respect  to  spatial  variables  as  well  as  time.  For  initial-boundary  value- 
problems,  however,  differentiation  with  respect  to  spatial  variables  “destroys”  the  bound¬ 
ary  conditions;  in  other  words,  we  do  not  know  anything  about  spatial  derivatives  of  the 
solution  on  the  boundary.  Since  the  boundary  conditions  are  important  in  the  integrations 
by  parts  leading  to  the  energy  estimates,  we  have  to  limit  ourselves  to  differentiating  with 
respect  to  time.  This  leads  to  energy  estimates  for  time  derivatives  of  the  solution.  Esti¬ 
mates  for  spatial  derivatives  must  be  obtained  in  a  different  fashion.  For  equations  such  as 
those  of  elasticity  this  is  done  by  using  elliptic  regularity  results  to  estimate  spatial  deriva¬ 
tives  in  terms  of  time  derivatives  (see  [3], [4], [8]).  In  the  present  paper,  we  show  how  this 
method  can  be  applied  to  the  Dirichlet  initial-boundary  value  problem  for  incompressible 
elasticity.  We  shall  prove  an  existence  theorem,  but  we  do  not  address  the  question  of 
whether  or  not  our  solution  is  the  limit  of  solutions  for  compressible  problems. 

We  consider  a  homogeneous  hyperelastic  body  with  reference  configuration  fl  and 
unit  reference  density.  By  x  we  denote  material  coordinates,  and  by  y(x,t)  we  denote  the 
spatial  position  at  time  t  of  the  particle  with  material  coordinates  x.  We  use  the  notation 
F'a  =  dy'/dxa  and  we  write  dxa  jdyx  for  the  components  of  the  inverse  matrix  F  1 .  The 
equation  of  motion  is  as  follows: 


^  dP  dx*  d*w(F)  ay 
V  dxa  dy 1  dF'adF30dxad x^  1  ‘ 


Here  W  =  W(F)  is  the  stored  energy  function,  p  =  p(x,t)  is  an  unknown  pressure,  and 
/  =  /(x,t)  is  a  prescribed  body  force.  Throughout  this  paper  we  adopt  the  Einstein 
summation  convention.  The  motion  must  satisfy  the  incompressibility  constraint 

det  F  —  l.  (1.2) 

We  seek  a  solution  to  (1.1)  and  (1.2)  subject  to  the  initial  conditions 

y(x,0)  =  y0(x),  y(x,0)  =  yi(x),  (1.3) 

and  the  Dirichlet  boundary  condition 

y(x,t )  —  x  for  t  >  0,  x  G  dfl.  (1.4) 

In  order  to  make  the  pressure  unique,  we  shall  normalize  it  by 


/  pOm) 

J  n 


dx  =  0. 


We  make  the  following  smoothness  assumptions: 

(Si)  fl  is  a  bounded  domain  (i.e.,  an  open,  connected  set)  in  R3  with  a  boundary  of  class' 
C3-1. 

1  We  could  weaken  this  assumption  and  impose  a  “Sobolev-type”  regularity  on  the 

boundary. 


(52)  The  stored  energy  function  W  is  of  class2  C5. 

(53)  The  initial  data  satisfy3  yo  €  H4{ fi),  yi  €  H3[Q). 

(54)  For  some  T  >  0,  we  have 


/€ 


f|  c*(|o,r);tf!-‘(n)), 


k= 0 


and  the  third  order  time  derivative  (in  the  sense  of  distributions)  of  /  lies  in 

L1((ofr];La(n)). 

Moreover,  we  assume  that  the  elastic  energy  satisfies  the  strong  ellipticity  condition. 


(E)  With 

A*0  =  d'W 
11  dF'adFf 


(1.6) 


we  have 

>  C|tl2|»P  (1.7) 

for  some  positive  constant  C  and  all  vectors  f,  77  6  R3. 

Finally,  we  have  to  assume  compatibility  of  the  initial  data  with  the  boundary  and 
incompressibility  conditions.  The  following  compatibility  conditions  are  required. 

(Cl)  y0  =  x  and  yi  =  0  on  dfl. 

(C2)  det  Vyo  =  1,  and  yi  satisfies  the  following  equation,  which  is  obtained  by  differenti¬ 
ating  the  incompressibility  condition  with  respect  to  time: 


dy\  dxa 
dxa  dy'Q 


(1.8) 


(C3)  The  initial  values  of  y  and  d3y/dt 3,  henceforth  denoted  by  y-i  and  y3,  vanish  on  dl 1. 

The  last  condition  requires  some  explanation.  The  initial  values  for  y  and  d3y/dt3 
have  to  be  determined  from  the  differential  equation.  However,  in  order  to  find  y  from 
the  equation,  we  must  first  find  the  pressure  at  time  t  =  0.  The  procedure  for  this  is  as 
follows.  First  we  differentiate  the  incompressibility  condition  once  more  with  respect  to 
time  to  obtain 

dy^dx^_dy^_dx^_dyJ_dx^_  = 
dxa  dyx  dxa  dyi  dx&  dy 1 


From  this,  the  quantity 


dy{  dxa 
dxa  dy1 


(1.10) 


2  We  always  think  of  W  as  being  defined  even  for  det  F  ^  1.  This  can  always  be 
achieved  by  extension. 

3  In  our  notation  for  Sobolev  spaces,  we  do  not  distinguish  the  number  of  components 

of  the  function,  i.e.  HA{ fl)  may  be  a  space  of  scalar-,  vector-  or  matrix-valued  functions 

depending  on  context. 


is  known  at  time  t  =  0.  If  we  now  apply  the  “divergence  operator”  (dx1  j dyx)(d / dxa)  to 
the  equation  of  motion  (1.1),  we  obtain  an  elliptic  differential  equation  for  p,  which  has 
the  form 

dx'1  d  rdxQ  dp  i 

dyx  dx'1  dyx  dxa  *  1’^1) 

where  g  is  known  at  time  t  =  0  and  lies  in  H1( Cl).  On  the  boundary  dCl,  we  multiply 
equation  (1.1)  by  (dx'1  /dyl)n^,  where  n  is  the  outer  unit  normal.  Since  (1.4)  implies  that 
y  must  vanish  on  3fl,  the  appropriate  boundary  condition  for  p  has  the  form 


dx*1  rdxa  dp  i  _ 
dyx  7  *•  dyx  dxa  ■* 


(1.12) 


where  h  is  known  at  time  t  =  0  and  lies  in  H3^2(dCl).  Equations  (1.11)  and  (1.12)  form 
a  Neumann-type  boundary  value  problem  for  p  (in  fact,  they  are  the  Neumann  problem 
if  we  transform  to  Eulerian  coordinates).  This  problem  has  a  unique  solution  p  6  H3((l) 
subject  to  the  constraint  (1.5),  provided  the  right  hand  sides  satisfy  the  condition 


/ 

J  30 


h  dS. 


(1.13) 


To  verify  that  (1.13)  holds  we  must  show  that  the  second  term  on  the  left-hand  side  of 
(1.9)  has  integral  zero  over  n.  If  we  transform  to  Eulerian  coordinates,  with  u  denoting 
to  the  velocity  field,  we  obtain  for  the  integral  in  question 


f  (Vyti)  :  (Vyu)T  dy  =  -  /  (u  -  Vj,)(divy  u)  dy  =  0. 

J  n  Ju 


(1.14) 


After  having  solved  for  p  at  time  t  =  0,  henceforth  denoted  p0,  we  can  insert  it  into  (l.l) 
and  compute  y  at  time  t  =  0.  Differentiating  (1.1)  once  with  respect  to  time,  we  obtain 
in  an  analogous  fashion  a  Neumann  problem  for  p,  which  we  must  solve  in  order  to  find 
d3y/dt3  at  t  =  0.  Condition  (C3)  means  that  the  so  computed  initial  data  for  y  and 
d3y/dt3  must  vanish  on  dCl. 

We  now  state  our  main  result. 

Theorem: 

Let  the  hypotheses  (Sl)-(S4),  (E)  and  (Cl)-(CS)  be  satisfied.  Then  there  is  a  T'  € 
(0,T]  such  that  the  problem  (l.l)-(l.5)  has  a  unique  solution  (y,p)  on  fl  x  0.Tr  with  the 
regularity  property 


4  4 

y  e  fl  c4 -*(jo.rj;//*(n)),  pef)  c4  *(  o.r  ,Hk  '(n)). 

k=0  *= 2 


4 


Remarks: 

1.  It  follows  from  the  Sobolev  imbedding  theorem  that  the  solution  satisfies  y  6  C2, 
p  6  C1  and  it  is  therefore  a  classical  solution  of  the  differential  equations. 

2.  The  result  can  be  modified  to  accommodate  Dirichlet  boundary  data  other  than  y  =  x 
provided  appropriate  smoothness  and  compatibility  conditions  are  satisfied. 

3.  The  result  can  be  extended  to  viscoelastic  materials  for  which  the  leading  order  terms 
in  the  differential  equation  are  like  those  in  the  elastic  case,  e.g.  K-BKZ  materials 
(see  [13]). 

4.  It  follows  from  (1.2)  and  (1.4)  that  the  mappings  x  — »  y(x,t)  are  actually  globally 
invertible,  see  [2]  and  [12]. 

5.  If  more  regularity  of  Q,  W  and  the  data  is  assumed,  and  the  appropriate  additional 
compatibility  conditions  are  satisfied,  then  the  solution  also  has  higher  regularity. 


2.  Outline  of  procedure 

The  solution  of  (1.1)-(1.5)  will  be  obtained  by  a  contraction  argument  which  involves 
solving  a  sequence  of  linear  initial-boundary  value  problems.  There  are,  of  course,  nu¬ 
merous  ways  to  formulate  (l.l)-(1.5)  as  a  fixed  point  problem;  for  technical  reasons  the 
precise  form  of  the  fixed  point  problem  is  a  very  delicate  issue.  Since  the  incompress¬ 
ibility  constraint  is  not  in  quasilinear  form,  we  shall  work  with  a  differentiated  version  of 
(1.1)-(1.5).  As  noted  in  the  introduction,  we  differentiate  with  respect  to  time  because  of 
the  boundary  conditions.  We  shall  actually  differentiate  (1.1)  twice  and  (1.2)  three  times. 
This  yields  a  linear  initial-boundary  value  problem  for  y  and  p.  Bounds  for  y,  y,  p  and  p 
will  be  obtained  via  elliptic  estimates. 

Assuming  that  (y,p)  is  a  sufficiently  smooth  solution  of  (1.1)-(1.5),  we  introduce  the 
notations  u  :=  y,  z  :=  y  —  A(y  —  x),  q  :=  p,  and  <f>  :=  p,  where  A  is  a  positive  constant  to 
be  chosen  later.  The  original  system  (1.1),  (1.2)  can  now  be  written  as  follows: 


a?  a*°  anr(Vy)  ay 

dxa  dyl  dFi'dFp  dxadx0 


(2.l)j 


det  (Vy)  =  1. 

After  differentiating  once  with  respect  to  time,  we  obtain 

dq  dxa  d2W{Vy)  <9V 


(2-1): 


z  —  — 


-  + 


dxa  dy'  dF^dFp  dxadx0 


d3w(Vy)  duk  ay  , 

dxa  dyk  dx0  dyx  dF^dF^dF*  dx'*  dxadxp  ^ 


du1  dxa 


dxa  dyx 
5 


=  0. 


(2.2), 

(2.2)2 


Differentiating  once  more  with  respect  to  time,  we  obtain 

d<t>  dxa  d2W (Vy)  d2z 


dx*  dy1  dx'1  dyk  dx “  dyx  dxa  dyk  ldx0 
d3W(Vy)  duk  d2u>  d3W(Vy)  ,dzk 


dF'a dF]pdFk  dx"1  dxadx0  dF'adF}pdFk  9x^  9x1 

a'wjvy)  du'  du'  ay 

dF^dF’pdF'dFl  dx~i  dx 6  dx°dx?  J 


dxadx0 


-  A**  -  A2y*  +  A  V  +  /*, 

(2-3), 

du‘  dza  du;  dz^ 

(2.3)2 

— - r  - -r - r  =  0. 

dz“  dyJ  dz^  dy‘ 

In  differentiating  terms  like  dxa /dy1  we  have  used  the  identity 

l-F"1  =  -F-^F-1.  (2.4) 

dt  dt  y  ’ 

Instead  of  using  (2.3)2  directly,  we  shall  differentiate  this  equation  once  more  with  respect 
to  time  and  apply  some  transformations  to  it.  In  doing  this,  we  use  the  identity 

i=W  =  0’  <25> 

which  is  obtained  as  follows.  Employing  (2.4)  with  d/dt  replaced  by  d/dxa,  we  find 

d  dxa  dxa  ay  dx P  j  dF.  dx^ 

dxa  dyl  dyi  dx°dx0  dyl  ^  dx0  dyl 


It  follows  from  (1.2)  that 


trF~'JS  =  0- 


Using  (2.5),  we  can  rewrite  the  first  two  terms  in  (2.3)2  as 

dx a  dy'  ;  dyl 

After  differentiating  with  respect  to  time,  this  expression  becomes 

d  \  dxa  dxa  ,  ■  dxa  duJ  dx0 

dx*  Z  Hy1  +  U  dy'  **  +  ^  ~  1  ^  dyl  dx 0  dy' 


*, ■ n\  ’T'> '/ V7/  V  V v  V  V  v  V  V  V  V  V  V  V  V  V  V  V  V  V*.*  V  *, 


For  dealing  with  the  last  term  in  (2.3)2,  it  is  convenient  to  transform  to  Eulerian  coordi¬ 
nates.  In  these  coordinates  the  term  is  simply  equal  to 


(Vj,u)  :  (Vvu)r  =  divj,  [(u  •  V„)u]  (2.10) 

(Vy  and  v  denote  Eulerian  space  and  time  derivatives,  while  ^  ’  +  (v  Vy)  denotes 

the  material  time  derivative).  In  (2.8),  we  set  v  =  (u  •  Vy)u  and  take  the  material  time 
derivative.  We  obtain 


$  +  (U  •  Vy)t>  -  (V  *  Vy)uJ  - 

To  proceed  further,  we  note  the  identity 

divy  [(u  •  Vy)v]  =  diVy  [«  diVy  t;  +  (v  •  Vy)u]. 


— divy  v  =  divy 


(2.11) 


(2.12) 


Using  this,  we  find 


v  =  divy 


t;  +  u  divy  v 


divy 


(u  •  Vy)u  +  (u  •  Vy)u  +  u  divy  v 


div. 


2(i  •  Vy)u  +  u  divy  v 


=  div. 


,  du 


2(^  ■  TV)U  -  2[((«  •  V,)«)  •  V„]»  +  u[(V,«)  :  (V,.)1] 


(2.13) 


By  transforming  this  back  to  Lagrangian  coordinates  and  putting  it  together  with  (2.9), 
we  fir. ally  obtain 


dy 


dxa  d{  r  ■  ,  du'  dx ^ 

-  —  lz'  +  Xu'  -  3  \z>  +  \(y>  -gt)]  —  — 
•  dxa  [  1  n  dx *  dy > 


kdx&  du1  dx'1  du'  ,  du 3  dx 0  duk  dx’1  j 
+  ^  dyk  dx &  dy 3  dx'1  U  dx P  dyk  dx ^  dy>  J 


(2.14) 


We  now  describe  the  basic  iteration  scheme.  Given  functions  y(n),  2(«),  P(n)  and  <?(n) 
on  Q  x  [0,X],  we  determine  y(„+i),  2(„+j),  P(n+i)  and  <^(„*i)  as  follows.  First,  we  set 
u (n)  :=  y(n)  and  :=  p(n).  For  technical  reasons  we  also  need  to  introduce  a  function 
U(n)  which  is  another  approximation  of  u  possessing  better  regularity  properties  than  U(n). 
We  set 


=  n^u(n),yi(i)  -  J  2(n)(x,  t)  +  A(y(n)(x,  r) 


(2.15) 


where  FI  is  the  projection  operator  from  an  appropriate  product  of  Hilbert  spaces  onto  the 
diagonal. 


We  then  determine  Z(n+i)  and  <£(n+i)  by  solving  a  pair  of  equations  related  to  (2.3)  i 
and  (2.14).  To  describe  this  pair  of  equations,  it  is  convenient  to  denote  the  expression  on 
the  right-hand  side  of  (2.3)  i  by 

G(y,  Vy,V2y,u,  Vu,  V2u,z,  Vz,/,x). 

We  replace  equation  (2.3)  i  with 

-  _ ^(n+i)  dx*  +  d2W(Vy(n))  d2z(„+p 
2(n+1)  dxa  dy'{n)  dF'adF]p  dxadx & 

+Gl{y(n),  Vy(n),V2y(n),u(n),Vu{n),V2u(n),z(n),V2(n),/,i).  (2.16)! 

In  the  analogue  of  (2.14)  we  shall  use  U(n)  rather  than  «(n)-  More  precisely,  we  replace 
(2.14)  with 

Jjfc  {*<-+'> +  Ai  w  - 3  W-o  + 

dlB  a“(»)  9x1  w;«)  .. 

dZ0  M»)  311  “<”>  dZ 0  3y(»)  ^  «Wn)  ' 

dxa  dr  1 

=:  =  0-  (2I6h 

If  we  impose  the  initial  conditions 


(2.16)2 


2(n+i)(x,0)  =  z0(x)  =  y2(x)-A(y0(x)-x),  i(n+i)(x,0)  -  x,(x)  =  y3(x)  -  Ay,(x),  (2.16)3 


the  boundary  cor  _n,ions 


and  the  normalization 


2(n+i)  =  0  for  x  €  dH,  t  >  0, 

/  <£(n+i)(x,0  dx  =  0, 

Jn 


(2.16). 


(2.16). 


the  problem  (2.16)  can  be  solved  uniquely  for  Z(n+i)  and  <£(n  +  1). 

Finally,  we  obtain  y(n  +  i)  and  P(n+i)  by  solving  a  nonlinear  elliptic  equation.  For  each 
t  t  j0,r),  we  consider  the  boundary  value  problem 

gjw,,  , 

3ia  dy!  anar)  r)Tadrp  vy(n^-i)  )  (r<-n  '  '  '  ' 1 


(n+l) 


dF^dF}p  dxadxp 


det  Vy(„.,)  =  1. 


(2.17)2 


(2.17)4 


y(n+i )(*,<)  =  X  for  x  e  dCl. 

Using  the  implicit  function  theorem  and  standard  results  for  linear  elliptic  equations,  we 
can  solve  this  boundary  value  problem  for  y(n+\)  and  p(n+,).  By  formally  differentiating 
(2.17)  with  respect  to  time,  we  can  obtain  estimates  for  y,  y,  p  and  p.  The  formal  differ¬ 
entiation  with  respect  to  time  can  be  justified  rigorously  by  taking  appropriate  difference 
quotients  and  passing  to  the  limit. 

To  implement  the  procedure  outlined  above,  we  must  solve  linear  initial-boundary 
value  problems  of  the  form 


2*  =  -X? 


dip 


+  A 


dV 


dxa  11  dxadx@ 


+  G\ 


(2.18), 


X, 


(2.18)2 


Ot0 

where  Gl,  Hl  and  the  coefficients  At]p  are  given.  The  solution  of  (2.18)  is  discussed  in 
section  4  and  will  be  based  on  the  Galerkin  method.  In  section  3,  we  discuss  the  needed 
results  for  elliptic  systems.  In  section  5,  we  use  the  contraction  mapping  theorem  to  show 
that  the  iteration  converges.  The  limit  will  provide  the  solution  we  seek. 

In  the  estimates  derived  in  the  following  two  sections  it  is  important  to  keep  track  of 
the  quantities  on  which  the  constants  appearing  in  the  estimates  depend.  In  the  remaining 
sections,  we  always  regard  ft,  W  (and  hence  the  constant  in  the  strong  ellipticity  condition), 
the  forcing  term  /  and  the  initial  data  as  given  once  and  for  all.  We  shall  therefore  not 
explicitly  point  out  any  dependence  on  these  quantities. 


3.  Elliptic  estimates 

In  the  proof  of  our  theorem,  we  shall  have  to  solve  linear  elliptic  systems  of  the  form 


dd>  dzJ 

(3.1), 

d!'  -  s 

X'  dl°  ' 

(3-1)2 

z  —  0  on  5ft, 

(3.1)3 

f  <t>  dx  ~  0. 

(3.1)4 

'n 


For  our  purposes,  it  is  particularly  important  to  understand  how  the  solution  (z,  o)  depends 
on  the  norms  of  the  coefficients  \ .  A°j \  B j'.  For  convenience  we  introduce  a  parameter 
K  which  controls  the  size  of  the  coefficients. 
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We  make  the  following  assumptions: 

(Si*)  The  coefficients  xf  and  are  in  Jf2(fi),  and  Bf “•  is  in  H 1  (fl).  Moreover,  we  have 
llxlli  +  M||*  +  ||B||,<if. 

(S2*)  For  some  given  integer  k  with  -1  <  k  <  1,  we  have  R  G  Hk(U),  S  G  iffc+1(fl). 

(I)  The  coefficients  xf  are  the  components  of  F~l,  where  F  is  the  gradient  of  a  globally 

invertible  mapping  y(x)  with  det  F  =  1. 

(E)  We  have  A°“f  =  A?°  and  the  strong  ellipticity  condition  (1.7)  holds. 

{C*)jugdx  =  0. 

The  following  lemma  holds: 

Lemma  1: 

Let  assumptions  (Si),  (Si*),  (S2*),  (I),  (E)  and  (C* )  hold.  If  A  >  0  is  chosen 
sufficiently  large  relative  to  K,  then  the  problem  (S.l)  has  a  unique  solution.  We  have 
z  G  Hk‘'‘(Q),  4>  G  Hk+x(U)  and  an  estimate  of  the  form 

IMU+2  +  \\<f>\\k+i  <  C{\\R\\k  +  ||S||*+I)  (3.2) 

holds.  The  constant  C  depends  solely  on  K  and  A. 

The  proof  employs  standard  techniques  in  elliptic  theory  and  we  omit  details.  One 
first  obtains  the  existence  of  a  unique  weak  solution  by  using  GArding’s  inequality  and  a 
standard  variational  argument.  Higher  regularity  is  obtained  by  using  the  fact  that  the 
system  (3.1)  is  elliptic  in  the  sense  of  Agmon,  Douglis  and  Nirenberg  [l].  The  regularity 
results  as  stated  in  [1]  only  apply  to  solutions  which  a  priori  satisfy  z  G  H2(Vl),  p  G  Hx((l) 
(as  stated  in  [l],  the  coefficients  would  also  be  required  to  be  of  class  C2,  but  it  is  clear 
from  the  proofs  in  [lj  that  this  is  not  actually  required,  since  H2((l)  is  an  algebra  in  three 
spatial  dimensions).  Regularity  of  weak  solutions  is  proved  by  Giaquinta  and  Modica  [6] 
in  the  case  of  traction  boundary  conditions.  As  they  remark,  their  techniques  can  also  be 
applied  for  Dirichlet  conditions. 

For  future  use,  we  note  that  if  /,  g  and  the  coefficients  (regarded  as  taking  values  in 
the  appropriate  Sobolev  spaces)  are  continuous  (or  bounded  measurable)  functions  of  a 
parameter  t,  then  this  property  is  inherited  by  the  solution. 

We  also  have  to  consider  the  solution  of  the  nonlinear  system  (2.1)  (with  the  boundary 
condition  (1.4)  and  the  normalization  (1.5))  for  given  z  and  /.  At  the  initial  time  <  =  0we 
have  the  solution  y  =  yo,  p  =  po  for  z  =  zo  (the  initial  value  of  z)  and  /  =  /o(x)  =  /(x,0). 
By  using  lemma  1,  the  implicit  function  theorem  and  the  regularity  theory  of  Agmon, 
Douglis  and  Nirenberg  [l],  we  obtain  the  following  result. 

Lemma  2: 

Let  A  >  0  be  chosen  sufficiently  large  (relative  to  the  data  for  the  original  problem).  If 
z  and  f  lie  in  // 1  (Q)  and  ||z-2o||i,  !|/~/o!!i  are  sufficiently  small,  then  (2.1),  (H),  (1.5) 
has  a  solution  (y,p)  G  H3(tl)  xH2(Cl).  Within  a  neighborhood  of  (y0,  po)  in  H3(U)\  H2(U) 
this  solution  is  unique  and  it  depends  smoothly  on  z  G  H  *(n)  and  f  G  Hl(Q).  If  moreover 
z  G  H2( D),  /  G  H2(U),  then  y  G  H4( H).  p  G  H3( H)  and  we  have  an  estimate  of  the  form 


v-TTTg  fv*VT.'r:  r*  r.-rj  wj  wjwjvj  wyrur'jwvwvwv  n  ws  tvw-jwvrju v  ~rjirj  <*-.  wj  irJ~T-"*\r%rjTrs  rj~v~j  rj 


where  the  constant  C  depends  only  on  A,  |j z  -  zo||i  and  || /  -  /o||i- 


4.  A  linearized  problem 


We  seek  a  solution  to  (2.18)  which  satisfies  the  boundary  condition 


2  =  0  on  dO, 

(4.1) 

initial  conditions 

2(z,0)  =  20(z),  i(x, 0)  =  2i(z), 

(4.2) 

and  the  normalization 

I  <fr(x,t)  dx  =  0. 

J  n 

(4.3) 

As  in  section  3,  it  is  imperative  to  understand  how  certain 

norms  of  the  solution 

depend  on  corresponding  norms  of  the  coefficients.  We  use  two  parameters  K  and  L  to 
measure  the  sizes  of  x?  and  A°f.  In  addition,  it  is  convenient  to  introduce  a  parameter 
U  which  measures  the  sizes  of  G  and  H . 

We  make  the  following  assumptions: 

(Si')  The  coefficients  x“  and  A°f  Wt  \n  W2<°°{\0,T},Hl{n))  n  W^flO.r];  2(0)) 
nL°°([0,T];  H3(Vl))  with  norms  bounded  by  L.  Their  norms  in  L°°([0,!r];  #2(0)) 
are  bounded  by  K. 

(S2')  G  e  C([0,T];L2(n)),  G  e  X#1  ([0,  T1];  X»2(0))  and  the  norms  are  bounded  by  U. 

(S3')  H  €  W2'1([0,T];L2(n))  n  I1([0,:T];  Jf2(n))  and  the  norms  are  bounded  by  U. 

(S4')  z0€H3(n),  Zi€H'( n). 

(E)  We  have  A°f  =  A?*  and  the  strong  ellipticity  condition  (1.7)  holds. 

(Cl')  zq  and  z\  vanish  on  dfl. 

(C2')  2]  satisfies  (2.18)2. 

(C3'j  H  vanishes  on  dO. 

(I)  The  coefficients  x“  are  the  components  of  F~l(t),  where  F{t)  is  the  gradient  of  a 
globally  invertible  mapping  y(x,t)  with  det  F  =  1. 

We  note  that  it  follows  from  (I)  that 


The  goal  of  this  section  is  the  following  lemma. 

Lemma  3: 

Let  assumptions  (Si),  (Sl')-(S4'),  (E),  (Cl')-(CS')  and  (I)  hold.  Then  the  initial- 
boundary  value  problem  defined  by  (2.18)  and  (4-l)~(4-S)  has  a  unique  solution  ( z,4> )  de¬ 
fined  onflx  [0,  TJ  with  the  regularity  property 

2 

2€  fl  Ck(\0,T]\H2~k{U)),  <peC(\0X\\H'(fl)). 


Moreover,  the  solution  obeys  the  a  priori  estimate 


E\z,<t>\(t)  <  T(U,K,T,T  •  L)  +  \(U,K,T,L)  [  E[z,0](s)  ds  Vt  €  \0,T],  (4.4), 

Jo 

and  hence 

E[z,<t>\(T)  <  T{V,K,T,T  ■  L)exp{T  A{U,K,T,L)),  (4.4)2 

where 

E\z,<t>\{t)  :=  max]^|||^(.,5)||^_*  +  ||^(.,a)||^.  (4.4)3 

r  and  A  are  functions  which  are  bounded  on  bounded  sets. 

It  is  important  to  note  that  T  depends  on  L  only  through  the  combination  T  ■  L,  and 
therefore  T  can  be  controlled  for  large  L  by  making  T  small. 

For  the  proof  of  the  lemma,  we  first  note  that  by  subtracting  from  z  the  reference 
function 

[  H(.,t)  dr 
Jo 

we  can  reduce  the  problem  to  the  case  H  =  0.  From  now  on  we  shall  therefore  assume 
H  —  0.  In  order  to  solve  (2.18),  we  differentiate  once  more  with  respect  to  time  and  then 
apply  a  Galerkin  method  to  the  time  differentiated  equation.  Differentiation  of  (2.18)  with 
respect  to  time  yields 

=  _Yaii  ,  A*i>  d2*j 

dt 3  dx«  dxadxP 


,a  *+  , 

X‘  dxa  '>  dxadxP 


+  G\ 


(4.5), 


adz'  .adz' 

XidP  +  X'd^  =  °‘ 


Because  of  (I)  the  latter  equation  can  be  rewritten  as 

+  *?*')  =  °  (4.S); 

Moreover,  recall  that  F £  is  the  matrix  inverse  to  x“-  We  can  write  (4.5) 2  in  the  form 

=°-  (4.5);- 

In  order  to  solve  (4.5),  we  must  prescribe  initial  data  for  z.  which  need  to  be  de¬ 
termined  from  (2.18) ,.  We  must  therefore  find  0  at  time  t  -  0.  For  this  purpose,  we 
apply  the  operator  xf  ^  to  (2.18),,  and  we  multiply  by  \fn^,  on  the  boundary  dfl.  In 
dealing  with  the  term  involving  z\  we  use  (4.5)2'  an<f  fact  that  z'  should  vanish  on  the 


boundary.  The  initial  values  of  z  and  z  are  known.  In  this  fashion,  we  obtain  a  problem 
of  the  following  form  for  <t>: 


=  X?  ; 


x?n**  x?j^  =  X?»V*‘  on  an. 


When  transformed  to  Eulerian  coordinates,  this  is  simply  the  Neumann  problem.  Subject 
to  the  constraint  (4.3),  there  is  a  unique  weak  solution  <j>  €  /f '(ft)  for  every  h  6  L2(ft)  (this 
is  known  as  the  Hodge  projection  theorem).  We  therefore  have  a  unique  initial  condition 
zi  G  L2(ft)  for  z  which  is  consistent  with  the  initial  data  (4.1)  and  boundary  data  (4.2). 
Instead  of  considering  (4.5)  \  directly,  we  shall  deal  with  the  equation 

d3z*  =  ,  A*fi  d2*’ 

dt3  *'  dxa  dx«dxV 

dd>  .  a  d3zi 

I"'7* 

where  <f>  and  z  are  expressed  in  terms  of  z  -  Az  by  the  following  elliptic  problem  obtained 
from  (2.18): 

+<a£|?-*5,  +  ci'  («)■ 

The  latter  equation  has  been  obtained  by  integrating  (2.18)2  (with  <7  =  0)  with  respect  to 
time.  Of  course  we  shall  eventually  show  that  z  =  z  for  the  solution  we  construct,  and 
hence  (4.7)  is  equivalent  to  (4.5)  j,  but  the  approximations  we  are  using  for  z  and  z  will 
actually  differ.  It  follows  from  the  elliptic  estimates  discussed  in  the  previous  section  that 
(4.8)  (with  the  boundary  condition  2  =  0  and  the  normalization  /n  4>(x)  dx  =  0)  can  be 
solved  for  2  and  <j>  if  A  is  chosen  appropriately.  Moreover,  an  estimate  of  the  form 

PIU~(|o,r|;ff*(n))  +  Plli-([o,r|;w*(n)) 

<  C^||z  -  Az||£,oo(|0,T];La(n))  +  il^ll  l°°  ( [0,T1;I-3  (o))  +  Po||//*(n)^  (4.9) 

holds. 

Let  now  V  be  the  space  of  all  divergence-free  vector  fields  in  Hq{CI)  and  let 
{£1^21  £3* •••}  he  a  basis  for  V.  Moreover,  let  V(f)  =  F(t)V  =  {u>  €  ffo(ft)  [F“'(f)u>  G 
V},  and  let  r?K(t)  =  F(t)£*.  Let  Vn  =  span  {£1,  £2, -.-,£»»}  and  Vn(t)  =  F(t)Vn.  We  seek 
an  approximate  solution  of  the  form 


*»( 0  —  Qn*(0>7<t(0, 


(4.10) 


(i.e.,  z  £  V„(t)),  which  satisfies  the  following  approximate  version  of  (4.7) 


33(Zn)’ 
dt 3 


Aapd2{ZnV 
t}  dxadz& 


."dUi «)» 

1  dxQ  dxadx & 


-  G*)  =  0 


(4.11) 


for  every  u>  €  Vn.  Here  4>n  and  zn  are  determined  from  (4.8)  with  z  on  the  left  side  of 
(4.8)i  replaced  by  zn ,  and  (.,.)  denotes  the  inner  product  in  L2( 0).  Finally,  we  define 
initial  data  by 


zn[xi^)  —  zno[x)  ~  zo[x)i  zn{xi  0)  —  Zn\(x)  —  P*Z\(x), 

5n(x,0)  =  zn2(x)  =  F(x,0)F_1(x,0)2ni(x)  +  P^22(x)-^(x,0)/'_1(x,0)2n,(x)^.  (4.12) 

Here  P,J  is  the  projection  in  V  (0)  onto  the  subspace  V„(0)  and  P°  is  the  projection  in 
L2(fl)  onto  Vn(0). 

Equation  (4.11)  is  a  linear  system  of  ordinary  integro-differential  equations  for  the 
coefficients  anK(t),  which  can  be  solved  by  standard  methods.  Hence  a  solution  to  (4.11) 
with  the  initial  data  (4.12)  always  exists.  To  obtain  an  energy  estimate  for  this  solution, 
we  set 

w  =  F~l  (zn  -  FF-'i^  =  ^  *„«(*)£«  (4.13) 

'  '  *=1 

and  integrate  from  0  to  t.  This  yields 

1  f  i=  /_  tu2  ,  ,afid(zny  d(zny 

iy0w*-‘»i  A->  ^~azTdz 

=  ^l^2lli»(n)  +  <f(<)fl(<)i»(<).^(0)-(F(O)f-'(O)2„„2n2> 

-  f  {F(t)F~'( r)in(r)  -  F(r)F-'(r)F(r)F-'(,)i„(r)  +  F(r)F->(r)J„(r),ln(r))* 

.'o 


l^^nl)*  »fid{zn  i)>  1  /’</a(2n)*  iQ^(in)^  J 

+  2(~d^’^  ~d^}  +  2  y0  -^r>  dT 


2  Jn  '  dx°  '  dx? 
9  A«0\d(zny 


-  f  dr 

J II 


+ 


dxa  %>  ’  dx & 

*P9{znY 


[  (  9  (ri  H(:  \k\  MaP9(zn)J  ,  . 

y0  'dx* 


+  /'«*»)•  -  +  Gl)  *.  (4.14) 

From  this  it  follows  that  we  have  uniform  bounds  on  the  norms  of  zn  in  L°°(l0,  T ];  L2(U)) 
and  zn  in  L°°(|0,  Tj;  H 1  (0)).  After  passing  to  a  subsequence,  we  may  therefore  assume 


weak**  convergence  in  those  spaces.  An  argument  along  the  lines  of  Chapter  III,  §1  in  [15] 
can  be  used  to  show  that  the  limit  z  is  in  fact  a  solution  of  (4.7)  and  (2.18)2.  Obviously, 
2  €  L°°([0,r];L2(n))  and  z  €  X>°°((0,  T’];  //" 1  (fl)).  Moreover,  it  can  be  shown  that  z  is 
weakly  continuous  into  L2( 0)  (and  in  that  sense  assumes  the  given  initial  data),  that  <j>  is 
in  L°°([0,r];L2(n)),  and  that 

d3z'  ,  a  dj> 


lies  in  I/°°([0,  T\;H~l  (fl)).  The  latter  need  not  be  true  for  each  of  the  two  terms  separately. 
The  reason  for  this  is  that  the  Hodge  decomposition  can  not  be  applied  in  H~x.  In  this 
respect  the  initial-boundary  value  problem  differs  from  the  problem  on  all  of  space. 

To  see  the  weak  continuity  of  2,  we  first  note  that  ( u>,FTz )  is  continuous  for  every 
solenoidal  test  function  u.  Therefore,  PFTz  is  weakly  continuous,  where  P  denotes  the 
Hodge  projection  in  £,2(fl).  Let  now  2  =  Fv.  Then  it  follows  from  (4.5)2  that 


v  -  F~lFF~1z  =  P(v  -  F-lFF~lz), 


-1  irr- 1 . 


(4.15) 


and  hence  v-Pv  is  continuous  in  L2(fi)  as  a  function  of  time.  As  a  consequence,  PFTFPv 
is  weakly  continuous,  and  since  the  operator  PFTFP  is  invertible  on  PL2(fl),  Pv  is  weakly 
continuous.  It  follows  that  v  and  hence  2  is  weakly  continuous. 

In  order  to  show  that  we  actually  obtain  a  solution  of  (2.18),  we  have  to  prove  that 
<t>  =  <t>  and  2  =  2.  By  integrating  (4.7)  with  respect  to  time,  we  find 


dz*  '>  dzadz0 


J0  dx°>  +  Aii' 


dxadx0  dxadx& • 


(4.16) 


In  view  of  (4.8),  we  can  rewrite  this  as 


d<t>  9^  .  a^,  dV  dzj 


_  _Z_)  +  AaP( 
Xt'dxa  dxa ’  A'}  ' 


dxadx0  dxadx0 


Jo  x'  dx°  dz°'  ,J  ‘ 


dzadx0  dxadx0 


)  dr  =  0. 


(4.17) 


Similarly,  by  integrating  (2.18)2  and  using  (4.8)2  we  obtain 


a/dz ‘  dz\  [*  .Q/dz*  dz\  J 

Xt  =  J0  *■  -  5?)  *■ 


(4.18) 


It  follows  from  the  elliptic  uniqueness  result  established  in  the  previous  section  and  a 
straightforward  perturbation  argument  that  z  =  z  and  <j>  =  <t>. 

It  remains  to  be  shown  that  2  is  actually  strongly  continuous  into  L2(U)  and  that  i 
is  strongly  continuous  into  Hl( fl).  To  do  this,  we  first  rewrite  (4.7)  in  a  different  form 


rc 


W 


SI 


by  first  multiplying  the  equation  by  FT  (this  transforms  the  <£-term  into  V<£)  and  then 
substituting  i  =  Fw  (this  transforms  the  incompressibility  condition  into  div  w  =  0.  We 
can  then  apply  the  technique  on  pp.  276-279  of  {llj,  which  employs  mollifiers  with  respect 
to  time.  This  is  why  it  is  important  to  write  the  incompressibility  constraint  in  a  form 
which  persists  under  such  mollification. 

Finally,  to  obtain  the  a  priori  estimate  (4.4),  we  proceed  along  the  lines  of  the  deriva¬ 
tion  of  (3.15)  in  [4].  The  basic  idea  is  that  by  virtue  of  the  weak-*  lower  semicontinuity 
of  L°°-type  norms,  (4.14)  yields  an  energy  inequality  for  z.  Using  this  inequality  together 
with  (4.9)  standard  Sobolev  estimates,  and  Gronwall’s  inequality,  we  obtain  (4.4).  See 
section  3  of  [4]  for  the  details  in  a  similar  situation.  Of  course,  in  the  derivation  of  (4.4), 
we  must  account  for  the  fact  that  we  have  subtracted  the  reference  function  H(.,t )  dr 
from  z. 


5.  Convergence  of  the  iteration 


In  this  section  we  prove  that  the  iteration  outlined  in  section  2  converges  in  an  ap- 

•  _  x  _  _  r  r  j*  rni  •.  .  •  •  it  i  i  i  e  r  *  •  /  ,  \ 


in  mis  section  we  prove  mat  me  iteration  outlined  in  section  i  converges  in  an  ap- 
propriate  space  of  functions.  The  iteration  involves  the  quadruplet  of  functions  (y,  z,p,<f>). 
Let  T  >  0  and  M  >  0  be  given  (later,  we  shall  choose  T  small  and  M  large).  We  denote  by 
\\y\\k,i  the  norm  of  y  in  W*,oo([0,  T]\ Jf^fi)).  By  Z(T,M)  we  denote  the  set  of  all  (y,z,p,^) 
which  satisfy  the  following  conditions: 


ML  +  llvlll,3  +  ML  < 

(5.1)i 

!IpIIo,3  +  l|p||?,s  +  IIpIII.i  <  M2, 

(5.1)2 

ii«»L  +  Mb  +  Mi. 0  <  m2. 

(5.1)3 

H^llo.1  < 

(5.1)4 

a(i,0)  =  z0{x), 

(5.2), 

y(z,0)  =  y0(i), 

(5.2)2 

det  Vy  =  1, 

(5-3) 

y  =  x  on  £H, 

(5-4), 

z  =  0  on  dU, 

(5.4)2 

u-t/e  L°°((o,rj; #3( n) n h^( n)) n w1'°°([o,:rj; h2( n) n h^{ n)) 

C  L2(|0,r];^3(n)  n  /f^(n))  n  //'(!o,r];//2(n)  n  h(\{ n))  =:  x,.  (5.5), 


.  v ,  -  -  -  • 
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On  the  other  hand  we  have 


yi(x)+  z(x,r)  +  \(y(x,r)  -  x)  dr 

Jo 

e  w'l‘“([o,ri;ff8(n)  n  H^(n))  nw2-°°(io,T];ff^(n)) 

cfrl((o,ri;/r2(n)n^(n))ni5r2(io,ri;fro1(n))  =.x2.  (5.5)2 

We  now  choose  II  in  (2.15)  to  be  the  orthogonal  projection  of  the  Hilbert  space  Xi  x  X2 
onto  the  diagonal  X\  n  X2.  Let  £  be  the  mapping  defined  by  the  procedure  of  section 
2  that  takes  (j/(n)t^(n)iP(n)i^(ti))  to  (y(n+i)i2(n+x)>P(n-n)»<^(n+i))>  It  follows  from  the 
results  of  section  3  and  4  that  £  is  well  defined  if  M  is  large  enough  relative  to  the  initial 
data,  T  is  small  enough  relative  to  M  and  A  is  sufficiently  large.  The  reason  why  T  must 
be  small  relative  to  M  is  to  ensure  that  ||zn+j(.,0  -  z0|li  remains  small  enough  so  that 
lemma  2  can  be  used  to  solve  (2.17). 

We  shall  show  below  that  if  M  is  increased  and  T  is  reduced  accordingly,  then 
(!/(»)>  2(n),P(n),<£(n))  converges  (in  an  appropriate  sense)  to  an  element  ( y,z,p,<j> )  of 
Z(T,M)  as  n  — ♦  oo.  We  shall  then  show  that  the  pair  (y,p)  is  a  solution  of  (1.1)-(1.5).  In 
the  proofs  of  lemmas  4  and  5  we  need  to  solve  the  nonlinear  elliptic  problem  (2.17)  as  well 
as  several  linear  elliptic  problems  of  the  form  (3.1).  The  parameter  A  must  be  large  enough 
so  that  these  problems  can  be  solved.  Examination  of  the  proofs  of  lemmas  4  and  5  reveals 
that  for  each  problem  of  the  form  (3.4),  ||x?j|2  +  +  ||2?,“||i  remains  bounded  as 

M  gets  large,  provided  T  is  small  enough  relative  to  M.  It  is  therefore  possible  to  select 
A  sufficiently  large  (relative  to  the  data  for  the  original  problem  (1.1)-(1.5))  so  that  no 
matter  how  large  M  is  made,  all  of  the  elliptic  problems  occurring  in  the  proofs  of  lemmas 
4  and  5  can  be  solved  provided  T  is  small.  We  now  fix  such  a  A  once  and  for  all.  Moreover, 
from  now  on  we  shall  always  assume  T  is  small  enough  (relative  to  M)  so  that  £  is  well 
defined. 

Lemma  4: 

If  M  is  chosen  sufficiently  large  and  T  sufficiently  small  relative  to  M,  then  £  maps 
the  set  Z(T,M)  into  itself. 


Proof: 

Let  M,T  >  0  and  ( y(n)»*(fi)iP(n),^(„ ))  given  and  set  (»(n+i),*(B+i),P(n+i)i 

</>(„■+ 1))  =  £(t/(n), 2(n), p(n), 0(n)).  Let  u(n),  U(„),  g(n),  G  and  H  be  as  in  section  2.  We 
first  apply  lemma  3  to  obtain  an  a  priori  estimate  for  2(n+i)  an<I  ^(n+i)-  For  this  purpose 

we  identify  A*f  with  ^j(Vp/n))  and  x“  with  -—f  .  Observe  that 

U^<\\£~^yo)h 


+ 


rT  PwjVyw)  a«f., 

o  dF^dFpdF*  dx 


(  ,0b  dt, 


(5.6) 
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from  which  we  conclude  that  ||A|!o,2  is  bounded  by  a  constant  plus  T  times  a  continuous 
function  of  M .  A  similar  comment  applies  to  ||xj*||o,2-  This  means  that  these  norms  can 
be  controlled  for  large  M  by  choosing  T  small. 

Next  we  note  that  the  appropriate  norms  of  G  and  H ,  i.e.  the  quantity  called  U  in 
lemma  3,  can  be  estimated  by  a  function  of  M  and  T  which  stays  bounded  if  T  is  chosen 
small  enough  relative  to  M.  The  reason  for  this  is  that  the  terms  appearing  in  G  and 
H  actually  have  better  temporal  regularity  than  is  required  by  lemma  3,  and  the  desired 
bound  can  be  obtained  using  Holder’s  inequality.  Employing  the  a  priori  estimate  (4.4), 
we  find  that 

E[z{n+1),<t>in+1)}(T)  <  r‘(M,r)exp(T.  A*(M,T)),  (5.7) 

where  T*(M,  T)  remains  bounded  if  T  is  chosen  small  enough  relative  to  M  and  A  *(M,T) 
is  bounded  for  M  and  T  in  bounded  sets.  It  follows  from  this  that  2(n+i)  and  <f>(n+i) 
satisfy  (5.1)3  and  (5.1)4  provided  M  is  large  and  T  is  small. 

To  obtain  estimates  for  y(n+i)  and  P(n+i)  we  proceed  as  outlined  in  section  2,  making 
use  of  the  results  of  section  3.  Observe  that 

II*(»+1)(mO-*o||i  <  /  \\z(n+i)(;T)hdT<MTVte\0,T}.  (5.8) 

Jo 

Thus  if  T  is  small  enough  relative  to  M  we  may  use  lemma  2  to  solve  (2.17)  for  j/(n+i)  and 
P(n+i)?  we  also  obtain  an  estimate  for  ||y(n+i)||4  +  ||P(n+i)||3  by  virtue  of  (3.3).  Formal 
differentiation  of  (2.17)  with  respect  to  time  yields  an  estimate  for  ||y(n+i)||3  +  ||P(n+i)||2 
by  virtue  of  lemma  1.  Of  course,  this  procedure  is  not  quite  legitimate  because  we  do  not 
know  a  priori  that  y(„+i)  and  P(n+i)  exist  as  functions.  However,  the  estimate  can  be 
obtained  rigorously  through  the  use  of  difference  quotients.  This  argument  can  be  applied 
once  more  to  obtain  a  bound  for  ||y(„+i)||2  +  IIP(n+l)||o- 
Putting  all  the  pieces  together,  we  conclude  that 

l|y(n+l)llo,4  *f"  lly(n-i-l)  II  1,3  +  II S/  (n-»- 1 )  II  2,2  +  l|P(n+ 1)  II  0,3 

+  ||P(n  +  l)lli,2  +  l|P(n+ 1 ) II  1,1  <  t(M,T),  (5.9) 

where  f  remains  bounded  as  M  gets  large,  provided  T  is  small  enough  relative  to  M .  The 
lemma  now  follows  from  (5.9)  and  our  previous  estimate  for  Z(n+i)  and  <£(n+i)-  We  note 
that  if  M  is  held  fixed  and  T  is  reduced,  then  £  will  still  map  Z{T,M)  into  itself. 

On  Z(T,M),  we  define  the  following  pseudometric: 

d{{y,z,p,<t>),(y,z,p,4>))  =  ||y  -  y||0,3  +  ||y  -  y||i,2  +  Up  —  p|io,2  +  Up  -  p|Ii,i 

+  H2  ~  ^11 0,1  +  li*  -  ^11 1,0- 

The  convergence  of  the  iteration  scheme  now  follows  from  the  following  lemma. 


(5.10) 


Lemma  5: 

If  M  is  chosen  sufficiently  large,  and  T  sufficiently  small  relative  to  M,  then  the 
mapping  £.*  Z(T,M)  — *  Z(T,  M)  is  a  contraction  with  respect  to  the  pseudometric  d. 


Proof: 

We  use  the  same  notation  as  in  the  proof  of  lemma  4.  In  addition,  vi  set  K(n)  := 
y(»+l)  -  !/(»)>  £(»)  ==  2(n+l)  “  2(n)i  P(n)  ■=  P(n+ 1)  ~  P(n ),  $(n)  ==  <f>(n+ 1)  “  <£(n),  G(n)  :  = 
G(y(n),Vy(n),V2y(R),u(n),Vu(n),V2u(n),2(n),V2(n),/,x),  H{n)  :=  if(u(n), Vu(n),y(n), 
Vy(n),2(n),  x).  We  now  fix  M  large  enough  so  that  lemma  4  applies  for  sufficiently  small 
T.  In  the  sequel  we  always  assume  T  is  small  enough  for  lemma  4  to  hold,  but  we  may 
have  to  choose  T  still  smaller.  Since  M  is  now  fixed,  we  suppress  the  dependence  of  all 
estimates  on  M.  It  is  convenient  to  define 

£;(!)  =  n»*(||Z(„)ll?  +  ll^(.)ll?+  lll'wlll  +  lln.)lll  +  IIJWI!  +  ll>(.)ll?)w-  (5.11) 
A  simple  calculation  gives  the  following  problem  for  Z(n+1). 


■7i  *♦<«■»)  .  any  .  *’*(.+■) 

(n+1)  <9x«  ayjll+1)  dFidF’g  l"+‘)  dxa3xf 


+  I 


r~{ 

(n+1)  V 


9<t>{n+l)  | 

f  dxa 

dxa  N 

dxa  1 

^(n+ 1)  ' 

d7W 

(V»(„)) 

\  a  <»+0 

dFtfFj, 

j  dxadxP 

fdZ\n+\)  d 

V  dxa  dxa 

(*(n+l)  - 

'(n+1)  _  G(n)>  (5.12)  i 


dyjn+1)V  dx°  dxa '  fn+»)  (n)/y 

/  dx*  dx •\(d*\n+l)  dHU)\_0 

V^y|n+1)  9y\n)J  \  dz°  dx<*  ) 


Z(n  +  l){xity  ~  ^(n+O^O)  — 

Z(n+i)(x,t)  =  0  for  x  e  dCl, 

/  $(n+1)(x,t)  dx  =  0. 

J  n 

We  now  use  (5.12)  to  obtain  an  estimate  for  ll^(»+i)lli  +  ll^(«+i)llo-  We  multiply  (5.12)  j  by 
Z(n+i)  and  integrate  with  respect  to  space  and  time,  performing  a  number  of  integrations 
by  parts  and  utilizing  the  remaining  equations  in  (5.12). 


(5.12) 2 

(5.12) 3 

(5.12) 4 

(5.12) 5 
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After  a  routine  (but  rather  long)  computation  we  find  that 

ll2(n+i)M)ll?  +  ll^(n+i)llo  <  r(r)£;(r)  +  a (t)  f  E'n+l(s)  ds,  (5.13) 

Jo 

where  T  and  A  are  continuous  functions  with  r(0)  =  0.  Differences  of  the  form  U (n)  := 
“(n+i)  -  «(n)  can  easily  be  estimated  in  terms  of  Y^n)  and  Z(n)  by  virtue  of  (2.15)  and  the 
fact  that  II  is  a  projection  operator.  The  reason  why  E’(T)  rather  than  E*  (t)  appears  on 
the  right  hand  side  of  (5.13)  is  the  nonlocal  nature  of  the  projection  operator  with  respect 
to  time. 

In  order  to  assist  the  reader,  we  give  the  details  of  a  typical  calculation  used  in  the 
derivation  of  (5.13).  In  particular,  we  shall  obtain  a  bound  for  the  integral 


■'U  \or^ar-p 

To  simplify  the  notation,  let  us  set 


(V»w)) 

la4.+.>, 

1 dxadx 

.,r))  dr. 

(5.14) 

'(v»(.))' 

){x,r). 

(5.15) 

We  first  apply  the  Cauchy-Schwarz  inequality  to  the  integrand  in  7(f).  Next,  we  observe 
that 

d2z ’ 

I!*?/  H°  5  c||*||J||*(.+1)||,  <  cm||*||,  (5.i6) 

by  standard  Sobolev  inequalities  and  the  fact  that  ||a(n+1)||2  <  M.  It  follows  from  the 
assumed  smoothness  of  W  and  standard  results  on  composite  mappings  that 


ll^llz  <CMl|VK(n)||2<cM||y(n)||3, 


(5.17) 


where  CM  is  a  constant  depending  on  M.  Employing  (5.16),  (5.17)  and  the  elementary 
algebraic  inequality  |afe|  <  A  (a2  4-  62),  we  find  that 

m\  <  \  j‘ WZ{.+  uUr)U  dr 

+  C  Jo  I|y(.)ll3  *  <c(tE;(T)  +  J‘e^(t)  dr)  Yf  £  [0.7-1.  (5.18) 

Since  M  has  already  been  fixed,  we  now  suppress  the  dependence  of  C  on  M  in  (5.18). 

To  obtain  a  bound  for  ||y(n+i)||3  +  !|P(n+j)jj2,  we  use  the  elliptic  problem 


dxa  dP 


(n+l) 


+ 


dy  in+!|  dx°  SF'adFi 


d"W  IV  »v 

-(vy(n+2))— ---  —  -  A Y, 


dx°dx 0 
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+ 


dFidF*'  **adzfi' 

det  Vy(n+2)  -  det  Vy(n+1)  =  1, 
yn+1(i,t)  =  0  for  x  6  dfi, 

I  P( n+l)  dx  =  0. 

Jn 


(5.19) 2 

(5.19) 3 

(5.19) 4 


The  terms  which  are  expressed  as  differences  in  (5.19)i  can  be  rewritten  in  the  form 
Bfj  f— r ,  using  the  identity 

v(Vv<„+,>)  -  V(V»,.+1))  =  (  j  vv(v»,„+1)  +  !jvy,„+I))  *,)  •  vy(n+1),  (5.20) 

which  holds  for  any  smooth  function  V .  A  similar  procedure  can  be  applied  to  (5.19)2- 
Applying  lemma  1  to  (5.19),  we  obtain  a  bound  for  ||y(n+i)ll3  +  ll^n+ijlh  in  terms  of 
||^(n+i)||i-  Similarly,  by  differentiating  (5.19)  with  respect  to  time,  we  obtain  a  bound  for 

11^(11+1)112  +  ll-P(n+l)||l- 

Combining  all  of  these  estimates  and  using  Gronwall’s  inequality,  we  find  that 


En+l{T)  <  r*(r)exp(r  •  A*(r))£;(r), 


(5.21) 


where  T*  and  A*  are  continuous  functions  with  T*(0)  =  0.  It  follows  from  this  that  £  is 
strictly  contractive  with  respect  to  d  if  T  is  sufficiently  small. 

From  lemma  5  it  follows  immediately  that  y(n),  Z(n)  and  p(n)  converge  in  the  metric 
given  by  d.  We  can  then  conclude  directly  from  (2.16)  1,  (2.16)5  that  <£(n)  converges  in  the 
sense  of  distributions  and  consequently  in  the  weak-*  topology  of  L°°[\0,  T);  H 1  (n)).  Let 
( y,z,p,4> )  be  the  limit  of  (y(n)»Z(n))P(n)><£(n)}  and  let  u  be  the  limit  of  U(n)-  It  >s  obvious 
that  (y,z,p,<(>)  satisfies  (2.16)  and  (2.17)  (with  the  indices  n  and  n  +  l  left  out).  To  verify 
that  we  have  in  fact  a  solution  of  the  original  problem  (1.1)-(1.5),  it  remains  to  be  checked 
that  y  =  z  +  A(y  —  x )  (this  implies  as  a  consequence  that  u  =  u)  and  p  —  4>.  For  this 
purpose,  we  first  note  that  by  differentiating  (2.17)  1  twice  with  respect  to  time  we  obtain 
(2.16)i  with  z  replaced  by  y  -  A(y  -  1)  and  4>  replaced  by  p  on  the  right  hand  side.  Next 
we  integrate  (2.16)2  with  respect  to  time.  This  yields 


dxa  dz' 
dyx  dxQ 


dxa  dz'0 
dyx  dxa 


dxQ  du3  dx &  dzx 


dy3  dxe  dyx  dxa  dyx  dxa 


f/*(u,Vu,y.Vy,z,x)  dr. 


(5.22) 
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By  differentiating  (2.17)2  twice  with  respect  to  time,  we  obtain  the  same  equation  with  z 
replaced  by  y  -  A(y  -  x)  and  u  replaced  by  u.  The  next  lemma  yields  an  estimate  for  u  -  u 
in  terms  of  z  -  y  +  A(y  -  x). 

Lemma  6: 

Let  X\,  X2  and  II  be  defined  as  in  (5.5)  above.  Then  there  is  a  constant  C  such  that 
for  every  u  €  X\  and  116X2  we  have 

||u  -  n(t»,t»)||Xi  +  || v  -  n(u,t;)||x3  <  C\\u  -  v||Hi([o,rj;H2(n)n//^(n))-  (5.23) 


For  the  proof,  we  note  that  it  follows  from  the  definition  of  n  that 
||u  -  n(u,t;)||^i  +  ||v  -  n(u,w)||^2  =  u_mm+J|r||^i  +  ||s|&3  =  ||u  -  vfXi+x^.  (5.24) 

Hence  we  want  to  show  that  the  norm  in  X1+X2  is  equivalent  to  that  in  H 1  ( [0,  T}\H2(Q)  n 
Ho  (ft)).  By  the  open  mapping  theorem,  it  is  sufficient  to  show  that  the  two  spaces  are 
equal  in  the  set  theoretic  sense,  i.e.  that  H  1(|0,r];H2(n)n//^(n))  is  contained  in  Xi  +  X2 
(the  reverse  inclusion  is  trivial). 

Let  A  be  the  Laplace  operator  on  ft  with  Dirichlet  boundary  conditions  and  A  =  —  A. 
We  have  /^(ft)  =  D{A1/2),  tf2(ft )  n  ff^ft)  =  D{A),  and  f/3(ft)  n  tf^ft)  D  D{A 3/2). 
Hence 

r,  =  L2(\0,T}]D(A3'2))n  H'(\0,T};D(A))  C  Xu  (5-25)! 

K2  =  H'(\0,T}-,D(A))nH2(\0,T}-,D(Al/2))  =  X2,  (5.25)2 

and  it  suffices  to  show  that  /f1([0,7'];  D(A))  C  Y1  +  Y2.  By  using  appropriate  extension 
operators,  we  can  further  reduce  this  question  to  the  analogous  problem  for  temporally 
periodic  functions  (with  any  given  period  larger  than  T )  rather  than  functions  on  the 
interval  |0,  T).  Hence  let 


Wi  =  L2p(\ 0, 2rr);  D{A3'2))  n  Hlp(\0, 2rr);  D(A)), 

W2  =  Hl(\0,2n\-D(A))nH2(\0,27r)-D(A^2)), 
W  =  Hp(\Q,  27r);  D(A)). 


(5.26) , 

(5.26)  2 
(5.26)  3 


The  subscript  p  indicates  that  we  are  dealing  with  periodic  functions. 

We  want  to  show  that  W  C  W,-t-H’2.  Let  <t>k  denote  the  normalized  eigenfunctions  of 
A  and  A|  the  corresponding  eigenvalues.  For  each  w  6  W  ,  we  have  an  expansion 


u-  =  2^  u’ki4>k(x)e,lt , 


IT;  w  N  (/  ~  -r  l)Aj  W 


(5.27) 


(5.28) 


■  .  .* 


.  '  .*•  .  .S,-' 


•  m  y  *  .  m  m  ■ 


*  ty.t 


The  desired  decomposition  is  now  obtained  by  setting 


,  l2  +  1  2  A?  , 

=  aj  +  /»  +  iw‘'-  *’*'  =  xl7W7Twl“-  (5  M) 

In  particular,  lemma  6  yields 

11“  ~  “  il  Ls  ([0,T];i/3  (fJ))  <  C|l!/  -  2  -  X(y  ~  x)  !liJ (|0,T);H2(n))-  (5.30) 

In  conjunction  with  the  uniqueness  statement  in  lemma  1  and  a  standard  perturbation 
argument,  this  implies  that  z  =  y  -  X (y  -  x)  and  <j>  =  p  as  claimed. 

Finally,  the  uniqueness  of  the  solution  of  (l.l)-(1.5)  follows  from  a  calculation  similar 
to  the  proof  of  lemma  5. 
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